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O Abstract. We present the results of theoretical and experimental studies of 

! l r dispersively coupled (or "membrane in the middle") optomechanical systems. We 

Ch calculate the linear optical properties of a high finesse cavity containing a thin dielectric 

membrane. We focus on the cavity's transmission, reflection, and finesse as a function 
of the membrane's position along the cavity axis and as a function of its optical loss. 
We compare these calculations with measurements and find excellent agreement in 
y—( cavities with empty-cavity finesses in the range 10 4 - 10 5 . The imaginary part of 

E> the membrane's index of refraction is found to be ~ 10~ 4 . We calculate the laser 

cooling performance of this system, with a particular focus on the less-intuitive regime 

in which photons "tunnel" through the membrane on a time scale comparable to the 

CO membrane's period of oscillation. Lastly, we present calculations of quantum non- 

demolition measurements of the membrane's phonon number in the low signal-to-noise 
<^ regime where the phonon lifetime is comparable to the QND readout time. 

OO 
O 



PACS numbers: 40.42, 42.50.-p, 42.50.Wk 



Dispersive optomechanics: a membrane inside a cavity 



2 



1. Introduction 



Nearly all the optomechanical systems which have been studied to date consist of an 
optical cavity whose detuning is proportional to the displacement of some mechanical 
degree of freedom. The mechanical degree of freedom is most commonly the position of 
the end mirror of a Fabry-Perot cavity [1-14] or the elongation of a waveguide [15-18]. 
In these systems the radiation pressure has a physically intuitive form: it is a force 
which acts on the mechanical degree of freedom and is proportional to the instantaneous 
intracavity optical power. 

Recently a new type of optomechanical system has been described in which the 
mechanical degree of freedom is a flexible, partially transparent object (such as a 
dielectric membrane) placed inside a Fabry-Perot cavity [19-24]. In this type of system 
the cavity detuning (and hence the radiation pressure) is periodic in the membrane 
displacement. 

Here we analyze several aspects of such a "dispersive" optomechanical device. We 
calculate its linear optical properties (transmission, reflection, and finesse) as a function 
of experimentally relevant parameters, and compare these calculations with experiments. 
We demonstrate a dispersive optomechanical device with a finesse F = 150, 000, and 
argue that it should be possible to realize F = 500, 000 using present-day technology. 

We also present calculations of the radiation pressure-induced cooling and heating in 
these systems. Because dispersive optomechanical systems consist of a compound optical 
cavity, their laser cooling is more complicated than in the more familiar "reflective" 
optomechanical devices described, e.g., in refs [1-18]. Lastly, we consider phonon 
quantum non- demolition (QND) measurements in dispersive optomechanical systems 
[19,25]. We focus in particular on phonon QND measurements with low signal-to- 
noise ratios, and consider how quantum effects might be manifest in such non-ideal 
experiments. 



2. Linear optical properties: calculations 

The geometry of the dispersive optomechanical devices considered in this paper is shown 
in figure [TJ Our one-dimensional model consists of two cavity end mirrors with electric 
field reflectivity r and transmission t (the two cavity mirrors are assumed identical 
in this paper, but the extension to unequal mirrors is straightforward). The dielectric 
membrane placed between the two end mirrors has a thickness and index of refraction 
n. The membrane's electric field reflectivity and transmission td are then given by [26] 

(n 2 — 1) sin knL c ] , , 

r d = !: L (1) 

2in cos knLd + {n 2 + 1) sin knLj 

■j- (^\ 

2in cos knLd + [n 2 + 1) sin knL^ 

where k is the wavenumber of the light incident on the membrane. 
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Figure 1. Schematic illustration of the dispersive optomechanical system. The 
membrane is represented as a mass on a pendulum. Enclosing the membrane are 
the end mirrors that define a high finesse cavity. Ai n , A re f, A tra n, A±, A%, A3, A4 are 
the incident, reflected, transmitted and circulating fields. 

Note that r d and t d are in general complex (reflecting the phase shift acquired 
by light reflected from or transmitted through a dielectric slab). If n is real, then 
I fd 1 2 + I td | 2 = 1. However in general n will be complex, with the imaginary part 
determining the membrane's optical absorption. 

To find the transmission and reflectivity of the cavity as a whole, we solve the 



following system of equations: 

A 1 = itA in + rA 2 e ikLl (3 a) 

A 2 = r d A ie lkLl + tt d A 4 e lkL2 (36) 

A 3 = it d A ie ikLl + r d A 4 e lkL * (3c) 

A A = rA 3 e lkL2 {3d) 

A refl =itA 2 e ikLl +rA in (3e) 

A tran = itA^ kL > (3f) 



Here A% through A4 are the electric field amplitudes of travelling waves in the cavity 
(as shown in figure 1), and A in , A re fi, and A trans are the amplitudes of the incident, 
reflected, and transmitted waves. L\ and L 2 are the lengths of the left- and right-hand 
halves of the cavity shown in Figure [T] [27]. 

Since we are primarily interested in cases where the cavity finesse is high and the 
membrane absorption is low, we find the cavity resonance frequencies by solving for the 
eigenfrequencies of the closed lossless cavity (i.e., assuming r = 1 and lm(n) = 0). The 
solution gives: 

4 0) = 20 r + 2 cos- 1 (| r d I cos 5) (4) 

where 8^' is the cavity's resonance frequency scaled by 2rc/ Jfsr an d Jfsr is the cavity's 
free spectral range. The scaled membrane position is 5 = 2k Ax. <p r is the complex phase 
of r d . From this expression it is clear that the magnitude of the membrane's reflectivity 
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Figure 2. Cavity detuning as a function of membrane position. The resonant 
frequencies of three of the cavity's longitudinal modes are plotted (in units of the cavity 
free spectral range) for several values of the membrane reflectivity. The membrane's 
power reflectivities | | 2 are the following: red: 0.000, orange: 0.080, yellow: 0.450, 
green: 0.773, blue: 0.982, and black: 0.999. For this calculation the various reflectivites 
were realized by fixing the membrane thickness and varying the index of refraction. 

| rd | determines the dependence of 5^ on the membrane position. <ft r , the complex 
phase of r d , sets an overall offset to 5^. 

The cavity detuning versus membrane position is shown for several membrane 
reflectivities in figure |2j For illustrative purposes the membrane reflectivity was varied 
by setting the membrane thickness to be unphysically small (Ld = 0.01 nm) and varying 
n. 

The analytic expressions for the transmission through the cavity and reflection from 
the cavity are also straightforward, but are too cumbersome to display here. We have 
not found a simple expression for the cavity finesse; instead we estimate it numerically 
from the linewidth of the transmission resonances. 

Figure [3] shows the finesse, resonant transmission (i.e., the transmission when the 
laser is resonant with the cavity), and the resonant reflection as a function of membrane 
position. These plots assume Ld = 50 nm, ^ = 1064 nm, r = 0.99991 and t = 5.28 x 10 -3 
(i.e., the power transmission of each end mirror is 16% of what it would be if it were 
lossless) corresponding to an empty-cavity finesse of 18, 000. In each plot the blue curve 
corresponds to n — 2.2 (i.e., a lossless membrane), while the red curve corresponds to 
n = 2.2 + 1.5 x 10 _4 z (i.e., a membrane with some optical loss). 

The differences between the curves for the lossless membrane (blue) and the lossy 
membrane (red) can be understood qualitatively. Placing a lossless membrane inside 
a cavity does not alter the rate at which photons leak out of the cavity; as a result 
the blue curve in figure [3|a) is flat. However the position of a lossless membrane does 
modulate the resonant transmission (figure |3](b)) because the cavity eigenmodes will be 
modified by the membrane. Cavity modes primarily localized on the right-hand side 
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Figure 3. Calculated finesse (a), resonant transmission (b), and resonant reflection 
(c) of the dispersive optomechanical cavity as a function of membrane position. The 
blue curves correspond to a lossless membrane, while the red curves correspond to a 
lossy membrane. Note that all the calculations assume the same (non-zero) loss in the 
end mirrors. This leads to a reflection and transmission which do not add to unity, 
and which depend differently upon the membrane position. The parameters used in 
this calculation are given in the text. 
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of the cavity in figure [T] will leak primarily out of the right-hand mirror, leading to an 
increased transmission coefficient for the cavity as a whole. Likewise, modes localized 
predominantly on the left-hand side of the cavity will couple primarily to external modes 
to the left of the cavity, leading to an increased reflection coefficient (figure [3]^c)). 

When n is complex, intracavity photons can be lost to the membrane absorption. 
This additional loss process lowers the cavity finesse by an amount proportional to the 
overlap of the cavity mode with the membrane, giving rise to the dips in the red curve 
of figure |3^a) when the membrane is positioned at an antinode of the cavity mode. Note 
that for the parameters used in this calculation, the finesse is not appreciably reduced 
from its empty-cavity value if the membrane is placed at a node of the optical field (i.e., 
corresponding to the peaks in the red curve of figure |3](a) . 

The loss of photons due to absorption in the membrane also prevents photons from 
transiting the cavity. As a result the resonant transmission has pronounced dips when 
the membrane is at an antinode (red curve in figure |3](b)). The reflection signal (figure 
|3]^c) arises from interference between intracavity light leaking out through the left-hand 
end mirror (which is affected by the membrane's absorption) and light promptly reflected 
from the left-hand end mirror (which is not). As a result its form is less intuitive, with 
membrane absorption leading to an asymmetric dependence on membrane position. 

3. Linear optical properties: measurements 

3.1. Experimental Setup 

Schematic illustrations of our experiment are shown in figures [4] and [5j Laser light is 
produced by a Nd:YAG laser (Innolight, Hannover, Germany) with wavelength A = 1064 
nm. The light passes through an acousto-optic modulator (AOM), and the first-order 
beam is sent to the optomechanical cavity via the steering mirrors Ml and M2. The 
cavity is formed by two dielectric mirrors each with a 5 cm radius of curvature (coated by 
Advanced Thin Films, Longmont, CO, USA). The mirrors are mounted to a cylindrical 
Invar spacer 6.7 cm long with a hole drilled along its axis to accomodate the cavity 
mode. 

The Invar spacer has a second hole drilled perpendicular to the cavity axis, allowing 
us to introduce the dielectric membrane into the waist of the cavity mode. The 
membrane is mounted on two piezoelectric elements. The first allows us to apply high 
frequency (~ 300 kHz), small amplitude (~ 1 nm) oscillations to the membrane (e.g., to 
excite its mechanical eigenmodes). The second piezo allows us to translate the membrane 
by roughly 2 /im along the cavity axis. The piezo elements are in turn mounted to a tilt 
stage (Thorlabs, Newton, NJ, USA KM05) which is rigidly attached to the Invar spacer. 
The tilt stage includes vacuum-compatible motorized actuators (Thorlabs, Newton, NJ, 
USA Z612V) allowing us to adjust in situ the angular alignment of the membrane 
relative to the cavity axis. In practice we have found that the membrane needs to be 
aligned to roughly 5 arcseconds in order to achieve the highest finesse described below. 
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Figure 4. Illustration of the optical cavity. The Invar cavity support (A) is mounted 
to the inside of the vacuum chamber. A series of alumina spheres (B) are mounted 
in cone-shaped recesses to provide kinematic mounting between the support and the 
Invar cavity spacer (C). The end mirrors (D) define the optical cavity. The membrane 
and piezoelectric elements are mounted to the Invar arm (E) which is in turn mounted 
to the tilt stage (F) . The tilt stage can be adjusted in situ by two motorized actuators 
(G). 

The membrane used in these experiments is a commercial, 50 nm thick, 1 mm x 1 
mm SiN x-ray window (Norcada, Edmonton, AB, Canada). The membrane is supported 
by a 200 /im thick Si frame. The exceptional mechanical properties of these membranes 
have been described elsewhere [28]. 

In practice we first align the cavity with the membrane removed. Then the 
membrane is inserted and its tilt and transverse position are adjusted until good 
transmission through the cavity is achieved. 

For most of the measurements presented here, the cavity is mounted inside a vacuum 
chamber which is pumped down to ~ 10~ 6 torr by an ion pump. Good vacuum is crucial 
to maintaining the membrane's high mechanical quality factor. 

3.2. Measurements 

We monitor the optical power reflected from and transmitted through the cavity using 
the photodiodes shown in figure [5j We can also determine the transverse profile of the 
cavity mode by imaging the transmitted beam with a video camera. 

Figure [6]^a) shows the optical power transmitted through the cavity as a function of 
the laser frequency and the membrane position when the laser is mode-matched to the 
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Figure 5. A schematic of the experimental setup. The membrane is depicted in grey 
between the two high finesse cavity end mirrors. The transmission and reflection are 
monitored with photodetectors (PD) and the signals are sent to a data acquisition 
(DAQ) card. 



cavity's TEM 0i o modes. The dark bands (indicating high transmission) correspond to 
the cavity's resonant frequencies. Comparison of these resonant bands with Q gives a 
value of | r c |= 0.35. As with all the data in this paper, the calibration of the membrane 
position is taken from the assumption that the features in the data are periodic in the 
membrane displacement with period rr/k. 

Figure |6](b) shows similar data, but taken with the input beam aligned in such a 
way as to couple to the TEM ,o mode and the nominally degenerate doublet consisting 
of the TEM 01 and TEM 10 modes [29]. In figure |6](c), the input is realigned to couple 
into still more of the cavity's modes. The cavity spectra shown in these figures can 
be easily explained using the standard description of higher-order transverse modes in 
optical cavities [30] , and a detailed description will be given in a future publication. 

The cavity finesse is determined from cavity ringdown measurements. In these, the 
laser frequency is swept slowly while the optical power transmitted through the cavity is 
monitored. When the transmitted signal exceeds a pre-determined threshold (indicating 
that the laser is coming into resonance with the cavity) the AOM switches off the input 
beam and the transient leakage of light out of the cavity is monitored. This decay has a 
single exponential form [19,28] whose time constant, r, is related to the cavity finesse, 
F, via F = 2nfFSRT- Figure [7|a) shows the finesse of the cavity's TEM ,o mode as a 
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Figure 6. Measured transmission (on a linear scale) as a function of laser frequency 
(in units of the cavity free spectral range) and the membrane position. In (a) the input 
beam is coupled almost exclusively to the cavity's TEM o (transverse Gaussian) mode. 
In (b) the input beam is coupled to both the TEMoo and the TEM m+n= i transverse 
modes. In (c) the input beam is coupled to several transverse modes, including the 
Gaussian. 

function of membrane position. The solid line in figure [7^ a) is a fit to the calculation 
described above and shown in figure [3]^a). 

Figures [7^b) and (c) show the cavity transmission and reflection on the TEM 0)U 
resonance as a function of the membrane position. The solid lines are fits to the data 
using the calculations described above. The data in Figure 7 were taken with the device 
in vacuum. 

The fits in figure [7] assume Ld = 50 nm, that the empty cavity finesse is 16,500, 
and that the end mirrors' transmission and reflection coefficients are t = 5.52 x 10~ 3 , 
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Figure 7. Measurements and fits of the cavity's finesse (a), resonant transmission 
(b), and resonant reflection (c) as a function of membrane position. The solid curves 
are fits to the data. These fits assume an empty-cavity finesse of 16,500 and give an 
imaginary part of the membrane's index of refraction of 1.5 x 10 -4 . 
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and r = 0.99991 (consistent with measurements of the cavity when the membrane was 
removed). Given these constraints, the fits yield n = 2.2 + 1.5 x lCT 4 z. The agreement 
between the data and fits indicates that our simple model does a reasonable job of 
describing the system. The few anomalous data points in figures (7^a), (b) and (c) 
correspond to membrane positions in which the TEM 0) o mode becomes degenerate with 
other cavity modes (see, e.g., figure [6](c)). Such degeneracies are not accounted for in 
our simple one-dimensional model. 

3.3. Discussion 

The loss in the membrane places limits on the maximum obtainable finesse. However the 
data in figure m(a) indicates that this limit depends strongly upon where the membrane 
is placed relative to the cavity nodes and antinodes. In order to determine this limit 
quantitatively, we use the value of Im(n) = (1.5 ±0.1) x 10" 4 extracted from the data 
and fits in figure |7[a) to calculate the finesse of hypothetical devices which are identical 
to the ones measured here but with higher reflectivity end mirrors. The result is shown 
in figure |8j and indicates that for experiments in which the membrane can be placed at 
a node, it should be possible to achieve F ~ 500, 000 with state-of-the-art end mirrors 
(i.e., those corresponding to an empty-cavity finesse of 1, 000, 000). For experiments in 
which the membrane must be placed away from a node, figure [7] (a) indicates that the 
optical loss in these membranes cannot be compensated for by better end mirrors. 

To test the prediction shown in figure [8] we replaced the end mirrors used in 
the measurements described above with end mirrors giving a measured empty-cavity 
finesse 205,000 ± 10,000. Figure [9] shows the resulting finesse (measured in air) as a 
function of membrane position, along with a fit which gives Im(n) = (2.3 ± 0.06) x 10 -4 
(the empty cavity finesse is set to 205,000). Although these fits indicate that there 
may be some sample-to-sample variation in the membranes' absorption, the overall 
level of agreement between the four sets of data (in figures [7] and [9J and the theory 
indicates that extrapolation to still higher-reflectivity end mirrors is justified, and that 
it should be possible to realize a dispersive optomechanical system with a cavity finesse 
~ 500,000 (see figure [8]). We have assumed L d = 50 nm throughout this discussion; 
still higher finesse could be achieved with thinner membranes, albeit with a decrease in 
optomechanical coupling. 

This high finesse is realized when the membrane is positioned at a node of the 
intracavity field. This arrangement is ideal for the phonon QND measurements described 
in [19] and below. However it also corresponds to a point at which the adiabatic radiation 
pressure is identically zero, and so may seem to preclude the realization of more familiar 
optomechanical effects such as laser cooling. In fact the situation is somewhat more 
complicated, as described in the following section. 
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Figure 8. A plot of the calculated finesse as a function of membrane position for 
different cavity end mirrors. All three curves assume the membrane's index of refraction 
is n = 2.15+1.5 x 10~ 4 i. The empty cavity finesse is taken to be 100,000 (lowest curve), 
314,000 (middle curve), and 1,000,000 (uppermost curve). When the membrane is 
positioned at a node, the device's finesse can be > 500, 000. 
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Figure 9. Finesse as a function of membrane position in a higher-finesse cavity. 
Similar to figure [3] (a), but taken with end mirrors giving an empty cavity finesse 
« 205, 000. The fit assumes an empty cavity finesse of 205,000 and gives Im(n) = 
2.3 x 10~ 4 . 



4. Optomechanical cooling and heating 

In this section, we obtain the optomechanical cooling and heating rates (as well as the 
optomechanical spring effect) in a membrane-in-the-middle (MIM) setup. We do so by 
solving the linearized coupled classical equations of motion for the membrane coordinate 
x and the amplitudes Ql,r of two optical modes in the left and the right halves of 
the cavity. This means we are assuming that the membrane reflectivity | r c | — > 1 
(corresponding, e.g. to the black curve in figure |2j and that the left and right half 
cavities are nearly degenerate (corresponding to membrane positions near the avoided 
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crossings in figure [2J). Although this regime has not been achieved experimentally, it 
is in this situation that the most striking deviations from the usual (i.e., "reflective") 
setup are to be expected and indeed are realized. Note that modeling a device in which 
| r c 1 (corresponding to an approximately sinusoidal relationship between 8j? and 
Ax) would require including the contributions from many different modes, not just two. 

To set the scene, we first review the calculation of the linearized optomechanical 
dynamics for the simpler, well-understood case of a single cavity mode in a "reflective" 
optomechanical device [31,32]. We will present it in a way that prepares us for the 
derivation involving the "dispersive" optomechanical device. 

4-1. Linearized dynamics of standard optomechanical systems 

We have the following equations of motion for the membrane coordinate x and the 
complex light amplitude a (rescaled such that a = 1 at resonance, and taken in a frame 
rotating at the laser frequency ujl): 

a = i(A — uj'x)a + — (1 — a) (5) 
x— — u>m(x — x ) — T M x + V\a\ 2 (6) 

Here A = ujl — uJ cauV [x = 0] is the frequency detuning of the incoming laser radiation with 
respect to the optical cavity mode frequency, to' = du cav /dx is the derivative with respect 
to the coordinate (u/ = —Ul/L in the usual setup, with L being the cavity length), k 
is the cavity's intensity ringdown-rate, ujm the membrane's mechanical frequency, xq its 
equilibrium position in the absence of light, and Ym its damping constant. The radiation 
pressure constant V introduced here has dimensions of frequency squared, and is given 
by V = —(oj'/u)L)E Tes /m , where m is the membrane's effective mass and E res is the 
light energy stored inside the cavity at resonance (proportional to the input intensity, 
E rcs = 4J in //t for a single-sided cavity). 

The steady-state solution (x,a) can be obtained by setting x(t) = x and a(t) = a 
and solving the resulting set of nonlinear equations. Note that for strong radiation 
pressure effects, more than one stable solution appears (two in the case discussed here). 
This is the static bistability that was found experimentally by Dorsel et al. [33] . We 
now assume this solution has been found and linearize the equations of motion around 
it, using \a\ 2 ss \a\ 2 + d*8a + c.c. 

8a = — iuj'(x8a + 8xa) + iA8a — —8a (7) 
8x = - u 2 M 8x - T M x + V(a*8a + c.c.) + f(t) (8) 

Here we have added a test force leading to an acceleration /. The response to this force 
will reveal the change in the membrane's damping rate and spring constant brought 
about by the radiation field. At a given driving frequency u , we decompose into 
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positive and negative frequency components: 8a(t) = 5a-e lwt + 5a + e +luJt , and likewise 
for Sx [where <5x_ = 5x* + due to 5x(t) being real- valued] . This leads to 

± iu5a± = —iuj'{x5a± + 5x±a) + iA8a± — — 5a±, (9) 

and therefore 5a± = Xa{^ zUJ )8x±, with the susceptibility Xa relating the light response 
to the membrane motion: 

= (A -u, +%)/<*>> -I (10) 

8b leads to 



— u 2 5x± = — uj 2 m 5x± =f iujT M 5x± + 

V(a*5a± + a5a* T ) + f± (11) 

After inserting 5a± and using 5x± = 5x^, we find the mechanical response 

5x ± = x(±w)/±, (12) 

where the mechanical susceptibility of the membrane has been modified due to the 
optomechanical coupling: 

X _1 (u;) = uo 2 M - u 2 + iuT M + S(a»). (13) 

All the novel effects are contained in the optomechanical "self-energy" 

E(co) = -V(a*xM + ax* a (-co)) . (14) 

The optomechanical damping rate may now be read off from the imaginary part of the 
susceptibility, evaluated at the membrane's resonance frequency: 

r opt = Im[E(u M ))/uj M , (15) 

which yields the known result for optical damping: 

r , V ._ |2 K J 1 1 1 ( . 

I op t = LO \a\ — < 2 2 ) 

2uJm 2 ( [oom - A + xuf + (f ) [-u M - A + xuf + (f ) J 
which is the difference between the rate of Stokes and anti-Stokes transitions. 

The prefactor is equal to — x\ V¥ u) 2 R nKj L 2 , where x| PF = h/(2mu)M)- Therefore, the 
optical damping rate is seen to obey the simple formula (see [31]) 

Topt = ^2 [Sff{wm) — Sff{— wm)\ j (17) 

where Spp is the spectrum of radiation pressure force fluctuations. 

The damping rate is positive at negative detuning (A — u'x < 0), corresponding 
to cooling, while it is negative at positive detuning, leading to an increase in the 
mechanical quality factor, parametric amplification, and, potentially, the onset of self- 
induced oscillations (once r opt + Tm < 0) [34]. 

Likewise, the shift of the mechanical resonance frequency (optical spring effect) is 
obtained from the real part: 

8u M = Re\E(u M )]/(2u M ) . (18) 
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4-2. Linearized dynamics of dispersively- coupled optomechanical systems 

We now turn to the "dispersive" optomechanical device. We model it by considering 
only two modes, residing to the left and to the right of the membrane. This is the 
correct description in the limit of a completely reflecting membrane. We consider the 
first deviation from that limit, i.e. the mode amplitudes cil and are coupled by 
photon tunneling through the membrane, at a frequency g (the tunneling amplitude). 
When the membrane moves to the right, the frequency of the right mode increases, 
while that of the left mode decreases. They are degenerate at x = 0, but the 
tunneling introduces a splitting and leads to new eigenmodes that are symmetric and 
antisymmetric combinations, as expected for any level anticrossing. 

These features are incorporated into the following equations of motion: 

r hl_ 

a = Ma+ 2 ? ( 19 ) 

where 



a = 


a L 


; M = 








have 


been 


rescaled 



i(A 



UJ X) 

-ig 



K-L 
2 



i(A + u'x) 



2 J 



(20) 



reach a value of 1 at resonance for illumination of the left/right cavity (though in 
the situation assumed here, the illumination is only from the left side, as displayed by 



the inhomogeneous term in 19)). A is the detuning of the laser from the (uncoupled) 
resonance at x — 0. We have assumed a real-valued tunnel coupling amplitude g. Note 
that the phase of an arbitrary complex amplitude g could be eliminated by incorporating 
these phases into the definition of cil and or. The optical resonances in the presence 
of coupling can be found by setting deiM(A) = (with kl = Kr = 0): 



(21) 



UJc 



^'cav,±W = ±\Jg 2 + (u'x) 2 

Comparing this with the general expression for the dispersion, 
(c/L) cos _1 (rdCos(47rx/A)), near the degeneracy point and for — > 1, we find the 
following relations to the dispersive device's parameters: 

g = (c/L)yj2(l-r d ) and J = -u L /(L/2), (22) 

where it should be noted that L is the full cavity length (comprising both halves), 



and y2(l — r^) ~ \td\ is the transmission amplitude of the membrane. Note that in 
expressions of this kind (like the one for to'), the optical frequency ujl is assumed to 
be that of the original modes at the degeneracy point (small deviations do not matter 
here). The membrane's equation of motion is of the form 



x 



-UJ 



M 



X 



Y M x + V{\a L \ 2 - \a R \ 



(23) 
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For simplicity, we have assumed the two halves of the cavity to be of the same length, 
which is the situation realized in the experiment (otherwise one would need to distinguish 
between uj' l and u' R as well as Vl and Vr). 

Once again, first the steady state solution is found from the system of equations 



a 







M' 1 



1 




u M (x - x ) - T M x + V(\a L 



l«fl| 2 ) 



(24) 
(25) 



Linearization around this solution and splitting into positive and negative frequency 
components as before leads to 

± iu)8a± = —iu'a z aSx± + A45a±, (26) 

where Ai contains x and a z is the Pauli matrix. Thus, we find 8a± = x a (±u)8x±, with 

Xa(u) = —uj'[iuJ — M\~ x io z oi . (27) 

The mechanical susceptibility, obtained by solving the linearized equation for 8x, is 
analogous to that found for the standard setup, see (13), except for containing radiation 
pressure terms from both the left and the right half-cavity: 



E( W )= -V(aix L M+&LXa(-u)) 

+ V(a* R x«(u J ) + aRX«*(-u)), (28) 

where x^J R refer to the two components of the vector Xa The optomechanical damping 



rate (and the spring effect) can now be obtained as before, from (15) and (18), by 



inserting (28). There is little point in displaying these lengthy expressions explicitly, so 



we will instead discuss the results in terms of plots, for the case of kl = kr = k. 
4-3. Discussion 

The diagram of damping rate r opt vs. membrane position x and detuning A (both 



measured in units of the optical resonance width) is shown in figure 10 r opt is 
determined by only two dimensionless parameters. These are the ratio of the membrane 
frequency Um to the cavity ringdown rate k and the ratio of the photon "tunneling" 
rate g to k: 

% , 9 
and — 

K K 

As long as the two optical resonances (i.e. the upper and the lower parts of the hyperbolic 
detuning curve) are separated by more than max(w M , k), they can essentially be treated 
individually. In that case, the behaviour of the damping rate r opt in the vicinity of each 
resonance is qualitatively the same as for a standard "reflective" setup. That means 
r opt is positive (negative) for laser light red-detuned (blue-detuned) with respect to the 
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Figure 10. (a) The cooling rate r opt for the membrane-in-the-middle setup, 
as a function of laser frequency detuning (with respect to the degeneracy point) 
and membrane displacement. Blue (red) corresponds to cooling (amplification). 
Parameters are g/n = 2 and Wm/k = 1- (b) Cross-sections of the preceding plot, taken 
for several membrane positions: x ■ \u>'\/k = 0, 0.5, 1, 1.5, 2, 4, 8. (c) Maximizing the 
cooling rate, as a function of membrane position, for various coupling strengths g. Note 
that the position has been rescaled by g, and the curves coincide for large g S> ojj\[, k. 
They saturate towards large displacements x, to the value given by the simple theory 
for the standard setup, while r opt vanishes at x = (see text) for large g. For small g, 
the maximum cooling rate can even become larger than the standard limit (see curve 
labeled g/n = 0.5). Note that r opt was maximized over the half-plane A < 0, where 
the global maximum of r opt is located when x > (the plot would be inverted with 
respect to x < if we were to maximize over A > 0). 



resonance, i.e. one has cooling (or amplification) for u>l < oj vcs {ojl > u> TeB ). When ujm/k 
is small, the maximum |r opt | is reached for a detuning of ±k/2 (the point of maximum 
slope in the intensity- vs-detuning curve). In the resolved-sideband regime uom 3> k, 
the maximum is reached at ±lum- The only quantitative difference is brought about 
by the change in the slope du res /dx, which is directly proportional to the net radiation 
pressure force acting on the membrane. As the slope goes to zero near the avoided 
crossing, so does r opt oc (duj ves /dx) 2 . Note that this is a result of our weak coupling 
approximation and linearization of the equations of motion. Inclusion of higher order 
terms would permit two-phonon Raman processes which can lead to cooling for red 
detunings of 2um- In addition, the circulating power is smaller when most of the light 
is stored in the right half of the cavity (since we assume illumination from the left), 
and therefore the cooling rate is correspondingly reduced on that branch (with positive 
slope du Ies /dx), as can be seen in figure 11 On the amplification side (r opt < 0), the 
membrane may settle into a state of self-sustained oscillations when r opt + T M becomes 
negative. Those regions of instability can therefore directly be read off diagrams such 
as those in figure 11 once the mechanical damping rate Vm is given. 

When the resonances touch, i.e. when they get closer than max(w M , k), the regions 
of cooling and amplification become visibly distorted, with intricate patterns as a result. 
In any case however, the diagrams remain inversion symmetric around the degeneracy 
point (upon simultaneous change of the sign of r opt ). 
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The physics of this regime can best be understood by analyzing the cases where 
the mechanical frequency Um becomes comparable to or even larger than the splitting 
2g of the dispersion relation, while the ring-down rate remains small. For the following 
discussion, we therefore refer the reader to the lower right panel of figure [TT| (g/x, = 1 
and uj m /k = 4). As seen in that figure, the cooling or heating rate is apparently 
maximal at places where the incoming radiation is in resonance either with the optical 
eigenfrequencies u> C£lVj ±(x), or with their sidebands u) csv> ±(x) ±ljm (full and dashed lines 
in that panel). The rate becomes particularly pronounced when these dispersion curves 
cross. At these places, there is interference between the eigenmode that is nearby in 
frequency, and the excitation of the other eigenmode via Raman scattering. Indeed, 
this interference is necessary to explain the remarkable fact that there can be some 
cooling or heating even at x = 0, particularly when 2g becomes smaller than u^: At 
the degeneracy point x = 0, we have the eigenmodes as symmetric and antisymmetric 
combinations of the basis modes, a± = (q^ ± an)/y2. Consequently, the radation 
pressure force 



Frad \oli\ — \ocr\ oc|a + + a_| — |a+ — a_| (29) 

vanishes identically, unless there is interference between ct + and a_. This can happen 
at the points where the dispersions u± )Cav (x) and lu± jCSCV (x) ±cjm cross. It also happens 
elsewhere, to a lesser extent, due to the finite cavity ring-down rates Kl,r, which broaden 
the resonances. The strongest effect is observed when ujm = 2g, where the resonance 



conditions are fulfilled simultaneously right at the degeneracy point (see Fig. 11, panel 



with g/ k = 2,uj m /k = 4; and Fig. 10[c), g/n = 0.5, uj m /k = 1). On the other hand, the 
effect vanishes in the limit ujm/9 0. In that limit, the standard picture is sufficient, 
when taking into account the suppression of the slope dcu cav (x) / dx at the degeneracy 
point, which reduces the cooling and heating rates to zero. 

In summary, in a certain regime it is possible to cool the membrane even at x = 
(i.e., where the bands are flat). This discussion may be important for experiments 
in which the membrane is kept at x = to realize a phonon QND measurement, as 
described in Ref. [19]. 



5. Signatures of quantum behaviour in a weak energy measurement 

As has been discussed above and elsewhere [19], a dispersive optomechanical device 
can be operated in a regime where the cavity frequency depends directly on x 2 , the 
position-squared of a macroscopic mechanical oscillator. As a result, one can make a 
direct measurement of the oscillator's energy E = Uu^i? 1 + 1/2) where n is the number 
of phonons in the membrane [19]. One drives the cavity on resonance and measures 
the phase of the transmitted beam; this phase is proportional to E. We thus have 
the possibility to detect a truly quantum aspect of the oscillator: the quantization of 
its energy. Note that this is impossible to do with a linear position detector (e.g. a 
cavity whose frequency depends directly on x), as in this case one measures both 
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Figure 11. The optomechanical cooling rate in the MIM setup near a degeneracy 
point. These plots show r opt as a function of membrane position \lo'\x/k (horiz.) and 
detuning A/re (vert.), for various values of lum/k and g/n. Blue refers to cooling 
(r op t > 0), red to heating or amplification (r opt < 0). 

the energy and phase of the oscillator, and is thus subject to the usual limitations 
imposed by quantum back-action [35]. All recent experiments in optomechanics and 
electromechanical systems have (to the best of our knowledge) employed linear position 
detectors, and are hence subject to these limitations. 

While having a non-linear coupling to the oscillator is certainly a prerequisite to 
detecting the quantum nature of its energy, it is not in itself enough: one also needs 
to consider the output noise of the detector. Here, this output noise consists of the 
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shot noise in the transmitted beam through the cavity, plus any additional technical 
noise associated with determining the phase of this beam. If this measurement was 
truly QND, this output noise would not be a problem: one could achieve any desired 
sensitivity by simply averaging the output signal for a sufficiently long time [35]. The 
back-action of a perfect QND measurement does not affect the measured observable, and 
thus the oscillator's energy would not fluctuate during the measurement. However, this 
is not the case for any real experiment, which is performed at some small but non-zero 
temperature. Because of this small temperature and the non-zero oscillator damping, 
the oscillator's energy will indeed fluctuate if one waits long enough. There is thus a 
limit to how long one can average, and thus to how well one can resolve the quantum 
nature of the oscillator's energy. 

In this section, we will address quantitatively the limitations on detecting quantum 
behaviour arising from the combination of the weak nature of the measurement and the 
unavoidable (thermal) fluctuations in the oscillator energy. We will focus on the realistic 
case where one can only obtain an energy resolution corresponding to a single quanta 
after averaging for a time comparable to (or longer than) the lifetime of a phonon 
Fock state. As such, one is no longer truly measuring the instantaneous energy of 
the oscillator; instead, one is measuring the time-integrated energy fluctuations of the 
oscillator. We will calculate this quantity for both a classical and a quantum dissipative 
oscillator, and will discuss whether differences between the two are experimentally 
resolvable. We will focus throughout on the experimental conditions proposed in 
Ref. [19]. In particular, we assume a situation where the oscillator is initially near its 
ground state, but is coupled to a dissipative bath with a temperature Tb a th hcuM/ks- 
In the proposed experiment this is realized by laser-cooling the membrane to its ground 
state and then shutting off the cooling beam while an energy measurement is made. 
Note that the somewhat related situation of QND measurement of qubit energy was 
studied theoretically in Ref. [36] . 

5.1. Model and measurement sensitivity in the zero-damping limit 

The quantity measured in the experiment is the phase shift of the transmitted beam 
through the cavity (or, equivalently, the error signal in a Pound-Drever-Hall setup); by 
dividing out a proportionality factor (the "gain" of the measurement), one can refer this 
signal back to the mechanical oscillator, expressing it as a number of quanta n(t): 



where n(t) is the actual number of oscillator quanta, and £(£) is the added noise of 
the measurement. We take £(£) to be Gaussian white noise with a (two-sided) spectral 
density S nn , i.e.: 



n(t) 



n(t)+m 



(30) 



{Z(t 1 )Z(t 2 ))=S nn 6(t 1 -t 2 ) 



(31) 
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For the cavity system and a shot-noise limited Pound-Drever-Hall measurement, one 
has: 

S - teXHl-rc) m 
* nn ~ 4096vrF 2 P in a4 1 ' 

Note that S nn depends both on the amount of output noise in the measurement, and 

on the strength of the cavity - oscillator coupling. 

We will be interested throughout in the case of a weak measurement, where one 

must time-average the output signal to counteract the effects of the added noise. We 

are thus led to the quantity m(t), the time-integral of the inferred number of quanta 

n{t): 



m(W) = / dt'n(t') = J dt'n(t') + J f (t') 

= + / m (33) 



As the added noise is white, the probability distribution P meas (m, ^avg) for the measured 
quantity m is a simple convolution of a Gaussian and the probability distribution P(m) 
of the true oscillator energy fluctuations m: 



-Pneas ^avg) / dm 

with 



P(m - m',t avg ) 



1 / (m'/t avg 

exp 1 



\2 N 



2VT ((T n t avg ) 2 V 2<J n 



(34) 




cr n (t avg ) = J - — (35) 



Consider first the ideal case where the measurement is completely back-action free, 
and where the oscillator damping 7 — > 0, meaning that there are no thermal energy 
fluctuations. In this case, if the oscillator starts with no quanta, it will always have uq 
quanta: P(m,t avg ) = 5{m — n t avg ). The distribution of m is then just a Gaussian: 

Pmeas(m, t avg ) = - = exp -5- (36) 

'27T {<J n t R ^ V 2(T ~ I 



To see evidence of the oscillator's quantum nature, we would like to be able to resolve 
Fock states that differ by a single quanta. For 7 = 0, these two states will each 
give Gaussian distributions of fh having means separated by t avg . As is standard, we 
can describe the distinguishability of these two Gaussians by a signal to noise ratio 
RsNR(ta. vg ). This is simply the ratio of the signal power to the noise power: 

d /, \ [(m)i - {m) 2 f C g ^avg ( „ 7 , 

K SNR{t avg ) —————— 2 - — ^2-7^ \ 67 ) 

(Ami + Am 2 ) {2t avg a n ) 4^ nn 

where Ami denotes the standard deviation of m for the first Gaussian distribution, etc. 

As expected, RsNR{t a vg) can be made arbitrarily large by increasing the averaging time 

t avg . In particular, the two Gaussians become resolvable (i.e. the averaged distribution 

has two as opposed to one maximum) when RsNR(ta, vg ) > ijf]. 

I Note that our definition of the SNR ratio is smaller by a factor of two than the SNR ratio S used in 
Ref. [19] 
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5.2. Distinguishing quantum from classical when 7 > 

The story becomes somewhat more complicated when we now include the unavoidable 
fluctuations of n(t). We will consider the experimentally relevant case where these 
fluctuations are only due to the dissipative bath coupled to the oscillator, and not to the 
back-action of the measurement. As discussed in Ref. [19], there is a small back-action 
effect associated with the fact that the cavity is coupled to x 2 and not the oscillator 
energy; this however is a much weaker effect than the thermal fluctuations we consider. 
The thermal bath coupled to the oscillator will cause a given oscillator Fock state |n) 
to decay at a rate T n . A simple golden rule calculation yields: 

r„ = 7 [n eq + n(2n eq + 1)] (38) 

where 7 is the damping rate of the oscillator, 

\ -1 



n eq = exp 



1 , (39) 



hu M 

A^Tbath 

and T bath is the bath temperature. 

Due to these thermal fluctuations, the distribution of the measured quantity rh will 
not be Gaussian. To obtain a very rough estimate of whether our measurement can still 
resolve quantum energy behaviour, we could still attempt to use the SNR ratio derived 
above; this was the approach taken in [19]. We assume that we start in the ground 
state (to maximize the Fock state lifetime 1/T n ), and use an averaging time equal to 
the lifetime of this state. Letting r = 1/T represent the lifetime of the ground state, 
we thus have as an approximate figure of merit : 

R = RsN R (t avg = t) = jf - = 1 (40) 

One might guess that if R > 1, one can resolve quantum aspects of the oscillator's energy; 
in [19], it was shown that achieving R ~ 1 could be possible in the next generation of 
experiments. However, the condition R > 1 is clearly an approximate one, as it neglects 
all the complexities arising from the fluctuations of the oscillator. In particular, the two 
distributions one is trying to distinguish are not Gaussian, and thus it is by no means 
clear that the SNR ratio R will remain a good measure of distinguishability. 

We will now assess more accurately the conditions required to resolve quantum- 
classical differences. The first step will be to ignore the added noise of the detector, 
and focus on the probability distribution of the "true" time-integrated oscillator energy 
m. We will do this in both the cases of a classical oscillator and a quantum oscillator; 
the respective distributions will be denoted P c i( m ;*avg) and -P q (m,t avg ). Having these 
distributions, we will then add the effects of the added noise S nn , and ask whether 
the corresponding measured distributions P me n S ,ci( m , ^av g ) and -P m eas,q( m 5 *avg) (as given 



by (34)) are distinguishable for a given level of noise and averaging time. As the 
distributions involved will be non-Gaussian, we will need to use a more sophisticated 
measure of distinguishability than the signal-to-noise ratio R used in the Gaussian case. 
We will make use of an information-theoretic measure, the accessible information X. 
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We start with the first step of our program: what are the probability distributions 
of the true integrated oscillator energy fluctuations ml Given the relative weakness of 
cavity-oscillator couplings, we will necessarily need to use averaging times comparable 
to or even longer than the lifetimes of oscillator Fock states. As a result, the experiment 
is no longer about measuring the instantaneous energy of the oscillator. Rather, we 
are asking whether one can see quantum behaviour in the energy fluctuations of the 
oscillator. The quantities we we wish to calculate (P q (m, i avg ) and P c i(m, i avg )) are thus 
formally analogous to the well-studied full counting statistics of charge in mesoscopic 
electron systems [37]; there, one wishes to calculate the statistics of the time-integrated 
current through a mesoscopic conductor. Given this similarity, we can employ a similar 
calculational technique in our problem. This was essentially done in Ref. [38], where 
the motivation was to describe an experiment where a qubit is used to detect Fock 
states in a nanoresonator. One calculates the dephasing of a qubit whose energy is 
directly proportional to the energy of a dissipative oscillator; this immediately yields 
the generating function of P q (m, £ avg ), P q (A,t avg ), defined by: 



P q (A,i avg )=/ dme~* Am P q (m,t avg ) (41) 

J — oo 

In Ref. [38], the focus was to understand the time-dependence of P q (A,t avg ), and hence 
the dephasing spectrum of the qubit. Here, the focus will instead be on its A dependence, 
as the Fourier transform of P q (A,t avg ) will yield the desired distribution P q (m, t avg ). 

Consider the initial condition corresponding to the proposed experiment: the 
oscillator is initially in a thermal state corresponding to a temperature Ti n ; t which differs 
from the bath temperature T bath . Using the method of Ref. [38], one finds that the 
corresponding generating function for a quantum oscillator is given by: 

P (A, t ) = e 7«/2 e -i(a-A)*/2 1 ~ M (42) 



with 



a = \J (X — ij) 2 — 4iA7n eq (43a) 

M = 2ArW ~ ~ A + ^ (436) 
2An ini t + (a + A - 27) 



Here, n eq is the Bose-Einstein factor associated with the bath temperature (cf. (39)), 
while riinit is the Bose-Einstein factor associated with the inital oscillator temperature 

T'mit- 

For comparison purposes, we also require P c x{m, t avg ), the distribution of m for 



a classical dissipative oscillator. To obtain this, we use -P q (A,t avg ) in (42) to find 
PE, q (s, t avg ) the distribution of integrated oscillator energy fluctuations s = Tiu{m + 
£ avg /2); this involves a simple change of variables. It is then straightforward to take 
the classical % — > limit to find Pe, c \( s , ^avg)- Defining the corresponding generating 
function via: 



iVOcW) = / dse-^ s P EA (s,t aYg ) (44) 
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we find: 



with 



P E ci(x,t) = e ^/2 e -^/2 1 M ci (45) 
' v ; 1 - M d e~ iat v y 



a c i = \l-l 2 - 4i/c B T ba thX7 ( 46 ) 
2k B T initX - (ad + nO 
2k B T init x + {a d - ij) 
The corresponding classical distribution of m, P c i(m, t), follows from: 

P c i(A, t) = P E ,d (x = ^ tj e M ' 2 (48) 

While the form of the classical generating function Pe,ci{Xi t) may seem unfamiliar, 
it is easy to check its behaviour in some simple limits. For example, consider the limit 
where 7 — > 0. Equation (44) then yields: 

PeAx) - . r 1 rr ( 49 ) 

which corresponds to a simple Boltzman distribution as expected: 

PeciM) = i—k — exp (- ^S 1 ) (50) 
' K ' k B T init t P V k B T init ) 1 ; 

We now have analytic expressions for the distribution of time-integrated energy 
fluctuations of both a classical oscillator (P c \(m,t), cf. (48)) and a quantum oscillator 
(P q (m, t), cf. (42)). The corresponding distributions of the measured quantity fa can 
be easily found by including the effects of the added noise S nn via (34); we denote these 
(respectively) as P mC as,ci(^) and P me as, q (m). To assess how different these two (non- 
Gaussian) distributions are from one another, we will consider their mutual information 
I. This measure of distinguishability is defined as [39]: 

T[P U P 2 ] = HHP, + P 2 )/2] - i (Pf [P x ] + H[P 2 \) (51) 

where H[P] is the Shannon entropy of the distribution P: 

H[P] := J dmP{m) log 2 P(m) (52) 

We can interpret the first term in Eq. (50) as the information in a signal in which 
each instance is drawn randomly from either Pi or P 2 . The second term is the 
average information under the same circumstances except that we are told from which 
distribution the signal is drawn. If telling us which distribution was used makes no 
difference then the two distributions are identical and the mutual information is zero. 
The larger the value of the mutual information X[Px,P 2 ], the more distinguishable are 
the two distributions Pi and P 2 . X is a convenient measure both because it is applicable 
to arbitrary distributions Pi and P 2 , and because it has a direct information-theoretic 
interpretation in terms of information transmission rates down noisy communication 
channels [39]. We note for two Gaussians distributions with identical standard 
deviations, X ~ 0.49 when Rsnr = 1, while X — > 1 for Rsnr. 3> 1- 
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Figure 12. Quantum jump traces in the presence of noise and temporal averaging, 
for different noise strengths ((a) and (b)) and increasing averaging time (top to 
bottom). The plots show the traces (n(t)) t that would be observed by doing a 
sliding time-average of n(t), including the noise £(t). Here the time-average was 
done by convoluting with a Gaussian whose width is set by t avg , in such a way as 
to have the simple relation = S nn /t avg . The noise strength has been chosen 

as S n n/T = 0.001 and S nn /r = 0.004 in (a) and (b), respectively. The time-interval 
displayed in each curve is 2r, where r is the ground state lifetime. From top to bottom, 
successive curves arise from the same trace by averaging over increasing time-intervals: 
t&vg/T — 0.01,0.05,0.1,0.2,0.5,1.0. Curves have been displaced horizontally and 
vertically for clarity. The values of a signal-to-noise ratio, defined in correspondence 
to the discussion in the main text, are (from top to bottom): (a) t& vs /(4:S nn ) = 
2.5, 12.5, 25, 50, 125, 250 and (b) 0.625, 3.125, 6.25, 12.5, 31.25, 62.5. 



5.3. Results 



Equations (34), (42) and (48) can now be used to quantitatively assess whether quantum 
versus classical differences can be resolved under a variety of different experimental 
conditions. In what follows, we will present only a few selected results relevant to the 
experiment proposed in Ref. [19]. 



5.3.1. Measurement runs starting in the oscillator ground state We first consider the 
ideal situation where the oscillator has been cooled to its ground state: T init = n init = 0. 
The cooling beam is then shut off, and the number-state measurement is made. During 
this time, the oscillator rapidly heats up due to its coupling to the equilibrium heat bath 
at temperature T bath . We will focus on the experimentally relevant case where 7 <C lom-, 
n eq 3> 1, and on averaging times small enough that the average number of quanta in 



Dispersive optomechanics: a membrane inside a cavity 



26 




50 100 150 200 50 100 150 200 

SNR Ratio (R) SNR Ratio (R) 



Figure 13. Left: distinguishability of quantum versus classical energy fluctuations, 
as measured by the mutual information X[P meaSiq , P meas c i], versus the inverse added 
noise R = l/(4S , „„r). In each case, we have used an optimal averaging time t avg which 
maximizes I and have started the oscillator in the ground state n init = 0. Following 
Rcf. [19], we have taken lu m /2it = 10 5 Hz, j/uj m = 1.2 x 10~ 7 and T bath = 300mK. 
Finally, we have shifted the quantum distribution so that both the quantum and 
classical distributions have the same mean, as in experiment, this shift of the mean 
would be hard to detect. Note that an SNR ratio of R = 1 in the Gaussian case 
corresponds to I ~ 0.49. Right: optimal averaging time versus R, same parameters. 



the oscillator remains much smaller than n eq . In this regime, the oscillator damping 7 
and the bath temperature Xb at h essentially only enter via the time-scale r, the lifetime 
of the n = Fock state. It is this time-scale which determines the initial heating-up of 
the oscillator: 

(n(t)> = n eq (l - e-*) (53) 
~ 77i eq t = t/r (54) 

Given these conditions, there are two relevant questions. First, given a certain noise 
level S nn , what is the optimal averaging time? Second, given that we have optimized 
the averaging time, how does the distinguishability depend on S nn l 

To illustrate this process, we show a series of calculated time traces corresponding 
to such a scenario in figure 12 for various values of the averaging time. In the left 
panel of figure 13 we show how the measurable distinguishability between the classical 
and quantum distributions depends on the noise level S nn , as parameterized by R 
(cf.(40)). The distinguishability is measured by the mutual information X between the 
expected experimental distributions for a classical and quantum oscillator (-P m eas,ci(^) 
and P m eas,q(^) respectively). Each run corresponds to starting the oscillator in the 
ground state and using an averaging time which maximizes X; the value of the averaging 
time is shown in the right panel of figure 13 Note that before computing the mutual 
information X, we have shifted the quantum distribution to remove the zero-point energy 
contribution (as resolving this difference in an experiment would be very difficult); the 
result is that both the classical and quantum distributions have identical means. We 
see that in general, one needs a noise level small enough that R ^> 1 to unambiguously 
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Figure 14. Distributions of the integrated output of the experiment, for both the 
cases of a classical oscillator (-P m cas,ci(?7i), dashed blue) and a quantum oscillator 
(P meaSj q(m), solid red) oscillator. Following Ref. [19], we have taken uim/Zk = 10 5 Hz, 
ll^M = 1-2 x 1CP 7 and Tbath = 300mK. In each case, the oscillator starts in the 
ground state, and an optimal averaging time has been used; we have also shifted the 
quantum distribution in each case to remove the zero-point shift in the average. The 
four panels correspond to different levels of added noise S nn : R = 200 (top left), 
R = 61 (top right), R = 11 (bottom left) and R — 1 (bottom right). In each plot, the 
range of m has been chosen to display 90% of the area of the quantum curve. 



resolve classical-quantum differences: the simple Gaussian estimate which suggests 
R ~ 1 is sufficient is too optimistic. Classical and quantum distributions of m for 



different values of R are shown in figure 14 



5.3.2. Starting at a finite temperature It is also interesting to ask what happens if the 
oscillator does not start in the ground state (i.e. n init > 0). In practice, one might not 
be able to cool the oscillator all the way to the ground state. Even if one could cool to 
the ground state, it would be very useful experimentally to extract as much information 
as possible in one run of the experiment. We saw above that if the oscillator starts 
in the ground state, the optimal averaging time is on the order of r. After this initial 
averaging time, the average number of quanta in the oscillator will be ~ 1. One could 
imagine starting a second averaging period at this point; the question is whether the 
initial temperature of the oscillator will make quantum versus classical differences even 
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Figure 15. Left: distinguishability of quantum versus classical energy fluctuations (as 
measured by T[P mca _ s , q , P moa s,ci]), versus initial oscillator temperature, and for a fixed 
noise strength corresponding to R = 11. Following [19], we have taken ujm/^tt = 10 5 Hz, 
j/ujm = 1.2 x 10~ 7 and Tbath = 300mK. In each case we have used an optimal 
averaging time, and have shifted each quantum distribution so that both the quantum 
and classical distributions have the same mean. Right: optimal averaging time as a 
function of ni n iti a i for the same choice of parameters. 



harder to see. 



In the left panel of figure 15 we consider a situation where the added noise 
corresponds to R = 11, and plot the distinguishability between P meas ,q and P mea s,ci 
(measured via X) as a function of the initial oscillator temperature n m \t- For each 
point, we have used an optimal averaging time; the dependence of this optimal time 
on ninit is shown in the right panel. As in previous plots, we have also shifted 
the quantum distributions so that they have the same means as the corresponding 
classical distributions. As could be expected, as the initial temperature increases, the 
distinguishability between classical and quantum distributions does indeed decrease. 
However, this decrease is slow enough that one could obtain useful information even 
if the oscillator starts at low but non-zero temperature. It is also interesting to note 
that while the overall distinguishability between the classical and quantum distributions 
decreases with initial temperature, the quantum distribution will develop multiple peaks; 



this is shown in figure 16 
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Figure 16. Distributions of the integrated output of the experiment, for both the 
cases of a classical oscillator (P m eas,ci(?7i), dashed blue) and a quantum oscillator 
(Pmeas,<i(m), solid red). Following Ref. [19], we have taken ujm/^tt = 10 5 Hz, 
ll^M — 1-2 x 1CP 7 and Tbath = 300mK. In each case, we have assumed a noise 
level corresponding to R = 11, used an optimal averaging time, and have removed the 
zero-point shift in the average of the quantum distributions. The panels correspond 
to initial oscillator temperatures of fimit — (top left), nj n i t — 1 (top right), n; n i t = 2 
(bottom left) and nj n it = 4 (bottom right). In each plot, the range of m has been chosen 
to display 90% of the area of the quantum curve. At higher temperatures, multiple 
peaks are visible, however the overall distinguishability from the classical distribution 
is reduced. 
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